In this study, we investigate the performance of the frozen-density embedding scheme within density-functional theory [ J. Phys. Chem. 97, 8050 (1993) ] to model the solvent effects on the electron-spin-resonance hyperfine coupling constants (hfcc’s) of the H₂NO molecule. The hfcc’s for this molecule depend critically on the out-of-plane bending angle of the NO bond from the molecular plane. Therefore, solvent effects can have an influence on both the electronic structure for a given configuration of solute and solvent molecules and on the probability for different solute (plus solvent) structures compared to the gas phase. For an accurate modeling of dynamic effects in solution, we employ the Car-Parrinello molecular-dynamics (CPMD) approach. A first-principles-based Monte Carlo scheme is used for the gas-phase simulation, in order to avoid problems in the thermal equilibration for this small molecule. Calculations of small H₂NO-water clusters show that microsolvation effects of water molecules due to hydrogen bonding can be reproduced by frozen-density embedding calculations. Even simple sum-of-molecular-densities approaches for the frozen density lead to good results. This allows us to include also bulk solvent effects by performing frozen-density calculations with many explicit water molecules for snapshots from the CPMD simulation. The electronic effect of the solvent at a given structure is reproduced by the frozen-density embedding. Dynamic structural effects in solution are found to be similar to the gas phase. But the small differences in the average structures still induce significant changes in the computed shifts due to the strong dependence of the hyperfine coupling constants on the out-of-plane bending angle.

The absorption spectra of aminocoumarin C151 in water and n-hexane solution are investigated by an explicit quantum chemical solvent model. We improved the efficiency of the frozen-density embedding scheme, as used in a former study on solvatochromism (J. Chem. Phys. 2005, 122, 094115) to describe very large solvent shells. The computer time used in this new implementation scales approximately linearly (with a low prefactor) with the number of solvent molecules. We test the ability of the frozen-density embedding to describe specific solvent effects due to hydrogen bonding for a small example system, as well as the convergence of the excitation energy with the number of solvent molecules considered in the solvation shell. Calculations with up to 500 water molecules (1500 atoms) in the solvent system are carried out. The absorption spectra are studied for C151 in aqueous or n-hexane solution for direct comparison with experimental data. To obtain snapshots of the dye molecule in solution, for which subsequent excitation energies are calculated, we use a classical molecular dynamics (MD) simulation with a force field adapted to first-principles calculations. In the calculation of solvatochromic shifts between solvents of different polarity, the vertical excitation energy obtained at the equilibrium structure of the isolated chromophore is sometimes taken as a guess for the excitation energy in a nonpolar solvent. Our results show that this is, in general, not an appropriate assumption. This is mainly due to the fact that the solute dynamics is neglected. The experimental shift between n-hexane and water as solvents is qualitatively reproduced, even by the simplest embedding approximation, and the results can be improved by a partial polarization of the frozen density. It is shown that the shift is mainly due to the electronic effect of the water molecules, and the structural effects are similar in n-hexane and water. By including water molecules, which might be directly involved in the excitation, in the embedded region, an agreement with experimental values within 0.05 eV is achieved.

We investigate the usefulness of a frozen-density embedding scheme within density-functional theory [ J. Phys. Chem. 97, 8050 (1993) ] for the calculation of solvatochromic shifts. The frozen-density calculations, particularly of excitation energies have two clear advantages over the standard supermolecule calculations: (i) calculations for much larger systems are feasible, since the time-consuming time-dependent density functional theory (TDDFT) part is carried out in a limited molecular orbital space, while the effect of the surroundings is still included at a quantum mechanical level. This allows a large number of solvent molecules to be included and thus affords both specific and nonspecific solvent effects to be modeled. (ii) Only excitations of the system of interest, i.e., the selected embedded system, are calculated. This allows an easy analysis and interpretation of the results. In TDDFT calculations, it avoids unphysical results introduced by spurious mixings with the artificially too low charge-transfer excitations which are an artifact of the adiabatic local-density approximation or generalized gradient approximation exchange-correlation kernels currently used. The performance of the frozen-density embedding method is tested for the well-studied solvatochromic properties of the n→π^* excitation of acetone. Further enhancement of the efficiency is studied by constructing approximate solvent densities, e.g., from a superposition of densities of individual solvent molecules. This is demonstrated for systems with up to 802 atoms. To obtain a realistic modeling of the absorption spectra of solvated molecules, including the effect of the solvent motions, we combine the embedding scheme with classical molecular dynamics (MD) and Car-Parrinello MD simulations to obtain snapshots of the solute and its solvent environment, for which then excitation energies are calculated. The frozen-density embedding yields estimated solvent shifts in the range of 0.20–0.26 eV, in good agreement with experimental values of between 0.19 and 0.21 eV.